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Abstract — Pulsed Radio Frequency Interference (RFI) has 
severely effected on communication system in common. Even 
though Low-Density-Parity- Check (LDPC) code which has the 
decoding performance to approach the Shannon Limit and it is 
designed as a powerful forward-error-correction (FEC) code in 
the Additive White Gaussian Noise (AWGN) channel, simulation 
results show that the performance of LDPC decoder is degraded 
when exposed to the pulsed RFI. According to such a pulsed RFI 
impact, joint iterative channel estimation and decoding technique 
is proposed in this paper so as to decrease the effect of pulsed 
RFI while less complicated in processing. The proposed method 
decreases the complexity by implementing the simple way of 
SNR estimation and improves the efficiency and applying joint 
iterative technique between channel estimation and decoding. 
For optimal decoder performance, Signal-to-Noise Ratio (SNR) 
estimation must be as accurate as possible. Because computed in 
each time of iterative decoder, noise variance has been always 
change resulting in the performance of LDPC decoder has better 
in efficiency. The simulation results show that both less complexity 
and higher outcomes can be achieved by this proposed approach. 

Keywords — LDPC decoding, Pulsed RFI, Pulsed RFI impact, 
Joint iterative, Channel estimation. 

I. Introduction 

Reliability is an important thing in communication systems. 
Undesirable signal should be eliminated. Not only planetary 
communications systems, but also satellite communications 
systems are commonly interrupted by pulsed RFI. Several 
researches which investigate the impact of pulsed RFI on 
LDPC decoder performance found that regardless of the level 
of SNR, the Bit-Error-Rate (BER) is always large [1]. Pulsed 
RFI frequently has the duty cycles in the order of 5% or 
smaller. This significantly decreases the impact on receiver 
performance compared to continuous interference with the 
same power and centre frequency [6]. In addition, pulsed 
RFI randomly turns up in the form of block and the noise 
variance varies with time. Currently, several applications, such 
as ISS ACS Transponders, 4th Gen Transponders, Integrated 
Receivers and satellite communication receiver, are involved 
with pulsed RFI, resulting in wide awareness of the impact of 
pulsed RFI. 

LDPC codes, invented by Gallager, achieve near capacity 
performance in a wide class of channels [5]. Generally, various 
simulations and applications are based on AWGN channel. 



However, our simulation is unable to consider being only 
AWGN channel because of pulsed RFI appearing. As men- 
tioned earlier, the noise variance of pulsed RFI is time varying 
so, SNR evaluation must perform on fly. 

Several researchers have investigated a number of SNR 
estimation algorithms [2]. Although all of them work very 
well when the SNR is high, their performance suffers at 
negative SNR. The iterative SNR estimation [3] is rest on an 
iterative solution for the maximum likelihood estimation of 
the amplitude from which the SNR is competent to compute. 
Its results show that it exhibits a lower bias and normalized 
MSE than other techniques and that the useful range extends 
to negative SNR. However, the iterative SNR estimation also 
has flaw because of demanding a lot of complicated processes 
and taking a large of calculating time depending on the number 
of iteration. In the co-channel interference scenario and non- 
AWGN channel, an initial Log-Likelihood-Ratio (LLR) has 
been developed in many ways [10-14]. This work constructs 
the initial LLR based on probability density function (pdf.) [4]. 
The different initial LLR for LDPC codes using the density 
evolution technique to analyze. Nowadays, joint iterative tech- 
nique have been extensively implement in several applications 
[15] [16]. The proposed method named joint iterative channel 
evaluation and decoding technique which has been developed 
to re-processing SNR estimation during each time of LDPC 
decoder processes in other words we combine working process 
between SNR estimation and LDPC decoder. 

The rest of this paper is organized as follows. Section 
II presents the system model, the pulsed RFI model and 
the influence of pulsed RFI on LDPC decoder performance. 
Section III, which reviews SNR estimation techniques, consists 
of second order moment SNR estimation and the iterative SNR 
estimation. In section IV, the joint iterative channel estimation 
and decoding approach is presented. The simulation results 
of the proposed technique which demonstrates in section V 
are compared to conventional and another complex decoding 
technique in term of BER. The last section is conclusion. 

II. System model and pulsed RFI model 

The whole system illustrates as Figure 1. and (1) where 
Xk is transmitted sequences and yk is received sequences. The 
transmitted signals (xk) are interrupted by pulsed RFI (n r fi^) 
and AWGN (n awgn ^) which both of them have zero mean, 
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but possess difference in noise variance. Distribution of AWGN 
and pulsed RFI are capable of representing as iV(0, a 2 = a 2 + 
Act) and iV(0, cr^) respectively. Practically, because most of 
the channel is time varying SNR mismatch over an AWGN 
channel, the noise variance is explained as o\ — a 2 + Act. 
When a 2 signify as the noise variance of AWGN and Aa is 
symbolic of the random walk signal of AWGN. 



Vk = Xk + n awgn ^ k + g(t) * n r f iik 



(1) 



where 

n aW gn,k ~ N(0,af ~- 

n r fi,k ~ ^(O^i) Figure 2. shows pulsed RFI in single 



+ Aa) 
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Figure 1. System model. 
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Figure 3. Pulsed RFI impact on LDPC decoder 
performance(the AR4JA (8192, 4096) code). 



source. Generally, pulsed RFI randomly turns up throughout 
the length of x k and has duty cycles 3-5% of the code length. 
In addition, pulsed RFI has the noise variance changing with 
time. Pulse train function (g(t)) is assumed to be periodic 

1 1 1 1 1 

0 r T 2T 3T 4T 

Figure 2. Pulsed RFI model. 



The iterative SNR estimation [3] is formed from the maximum- 
likelihood principle that has even lower estimation bias at low 
and negative SNR. They develop an iterative search algorithm 
to find the amplitude that satisfies (4) by defining function 
F(x) as (5). They must perform the iterative algorithm to find 
the root of F{x) where F{x) = 0. the SNR is able to estimate 



as 



1 N 



Ayk, 



(4) 



k=l 



where th(x) = (e x — e x )/(e Xj ce x ) 



gating function with pulse repetition rate, rate 1/T and duty 
cycle r/T [1]. Figure 3. demonstrates the influence of pulsed 
RFI on LDPC decoder performance. It is found that whatever 
SNR is, BER under pulsed RFI environment is always high. 

III. SlGNAL-TO-NOISE RATIO (SNR) ESTIMATION 
TECHNIQUES 

Various digital communication applications, such as power 
control, bit error estimation, and turbo decoding, involve 
the knowledge of the SNR. For optimal performance, SNR 
estimation must be as accurate as possible. Several techniques 
have been proposed for SNR estimation. In this work, two 
techniques which are taken into the simulation. 

The second order moment SNR estimation [2] is popular 
in several applications because it provides a simple processing. 
SNR is capable of computing from the ratio of the signal mean 
(A) squared and the noise variance (a 2 ) which both of them 
are expressed as 

1 N 

^kl (2) 



k=i 

1 N 



N 



Af 



(3) 



F{x) 



1 N 



xyk , 

st2 1 



(5) 



k=l 



Therefore, the complexity in computation of two methods have 
significant difference. In other words the second order moment 
SNR estimation technique is less computational processing 
than the iterative one. 



IV. Joint iterative channel estimation and 

DECODING 

The system model which mentioned in the section II points 
out that both AWGN and pulsed RFI are capable of explaining 
as the normal distribution. Based on the pdf of the distribution, 
we are able to form the initial LLR of simulated channel. If 
the assumed noise pdf has distribution as 



/(«) 



where <j\ > a 2 



y^ 2 



+ 



(6) 



k=l 



Contaminated Gaussian Log-Likelihood-Ratio (CGLLR) 
[4] has been constructed. CGLLR provides initial information 
more accurate than implementing Gaussian LLR. CGLLR is 
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re-written as the follow 
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Where 

• a\ = a 2 + Act, is the noise variance of Gaussian part 
of the pdf. 

• a 2 is the noise variance of the contaminating heavy- 
tailed pdf. Chosen ranges are 1.0 < o\ < 4.0. 

• e is the percentage of sample from the heavy-tailed 
pdf. e can range from 0 < e < 0.5 with value greater 
than 0.3 occurring with very low probability. 

To implement this kind of LLR equation, two parameters must 
be computed. 

The first, the percentage of sample from the heavy-tailed 
pdf (e) or pulsed RFI part. Algorithms that are based on higher 
order moments are used to detect the presence of pulsed RFI 
[7] [8]. Figure 4. illustrates detection of pulsed RFI technique 
by second order moment. Let N is code length and G is 

pulsed RFI area 
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Figure 4. Detection pulsed RFI by the second order moment. 



number of group. N is divided into G groups; so, each group is 
composed of N/G samples. Because the second order moment 
(M 2 ) value of pulsed RFI area is higher than any other region, 
we take the advantage from this point to capture the area that 
pulsed RFI appearing. The second order moment (M 2 ) of every 
group is calculated and the group which has the largest M 2 
value is capable of considering to be pulsed RFI region. Hence, 
number of sample of pulsed RFI is approximately equal to 
N/G. To ensure that all pulsed RFI samples are detected, a few 
samples both left and right side of that group are included. In 
other words the pulsed RFI region is approximately equal to 
the area of group which it has the highest M 2 value including 
its side-band. Further precision may divide N into several 
groups and each group contains smaller samples. Applying 
the same technique; find M 2 of every group then pulsed RFI 
area may contain more than one group which have the highest 
M 2 value respectively. In practice, the second order moments 
(M 2 ) are estimated by their respective time averages for both 
real and complex channels as 



M 2 



1 N 



N 



(8) 



k=i 



The second factor that needs to compute is the noise variance 
of Gaussian part of the pdf (a 2 = a 2 + Act) because the 
channel status information is unknown. In various applications, 
the second order moment SNR estimation is used because it 
is the simple way to obtain o\ and does not take time to 



calculate a lot. However, the channel in this work is not same 
as other works since pulsed RFI also appears in the channel. 
In this environment, the second order moment SNR estimation 
is unable to provide very accurate value. The iterative SNR 
technique offers the noise variance value very near to the 
ideal value. Therefore, LDPC decoder performance has an 
exceptional efficiency. However, the iterative SNR estimation 
demands a large amount of processing time for the result 
depending on the number of iteration. 

To obtain higher performance than conventional method 
which deployed the second order moment SNR estimation and 
less processing time than the complex method which imple- 
mented the iterative SNR estimation, the proposed method has 
been created. 

Let consider in case of all BPSK modulated information 
in real channel is 1 or -1, the probability density function of 
a received vector (yi, y 2 , . . . Vn) is expressed as 



N 



fn(yi,y2,:-,yN) = Yl 



l 



27TO- 



(9) 



where transmitted signal taking values from{A, —A} and let. 

dfn(yi,y2,-..,yN) 



Therefore, 



dA 

dln[f n (y u y 2 ,...,yN)] 
OA 



= 0 



0 



(10) 



(11) 



From (9) and (11), the noise variance is computed 

N N 

- 2 = ^£^-^£^) 2 w 



Giving concentration on (2), (3) and (12), it is found that both 
the signal mean value (A) and the noise variance (a 2 ) of this 
case are as same as the second order moment SNR estimation. 
Therefore, the second order moment SNR estimation technique 
is the best solution to evaluate SNR of this kind modulated 
information. 

Even though initial information are randomly generated, 
we are able to modify modulated information to be the new 
information which contains only 1 or -1 as same as the former 
case. The proposed approach is constructed from this principle. 
Joint iterative channel estimation and decoding are generated 

Iterative Times 
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Pulsed RFI 



SNR 
Estimation 



-LcG(Yk)- 



LDPC 
Decoder 



L CN&VN update- 



Figure 5. Schematic block diagram of the proposed 
technique. 
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based on the fundamental that each iterative time of decoding 
process, the noise variance is re-computed by the second 
order moment SNR estimation of the modified information. 
Schematic of the proposed technique shows as Figure 5. 

Taking the advantage from the second order moment (M2) 
as mentioned earlier, the received signal (yk) is erased pulsed 
RFI area. The received signal after deleting pulsed RFI region 
denote as yk_noRFi- The second step is the noise variance 
(af) evaluation. On the one hand, if it is the first time of 
LDPC decoder iterations, yk_noRFi will be used to perform 
the noise variance calculation. On the other hand, the noise 
variance is computed from yk_new For optimal LLR, LcG(Vk) 
is determined in (7). We carry out the LDPC decoder only once 
iterative time. The results consist of of hard-decision output, 
values of check node (CN) and variable node (VN) updating. 
Hard-decision outputs are feed back to multiply with the pre- 
eliminated pulsed RFI signal (yk) resulting in the new received 
signal (yk_new)- In order to keep LDPC decoder processes, 
values of check nodes (CN) and variable nodes (VN) updating 
are also feed back so as to use in the initial state of the 
next LDPC iteration. In the check stopping criteria step, if the 
number of iterative time is not equal to the maximum LDPC 
iterations or decoding error is not equal to zero, go back to 
step 2.) and increases the number of iterations by one until 
reach the stopping criterion. 



procedure proposed approach (y k ) 
Eliminate pulsed RFI. 

if decoding iter. — 1 then 



4: 



else 



var(y k _noRFi) 
var(y k _new) 



5: end if 

6: LLR computation 

7: Once iterative time LDPC decoder 

8: Output feedback 

Vk_new = Vk ' hard - decision 
C heck N ode (C N) updating value 
V ariableN ode(V N)updating value 
9: Check stopping criteria : 

10: if Iter, times = max. iteration or Error — 0 then 
Go to 13. 

11: else 

Back to 3. 

12: end if 

13: Output codeword 

14: end procedure 



V. Simulation results and discussions 

The LDPC code studied here is the AR4JA (8192, 4096) 
code recommended by the Consultative Committee for Space 
Data Systems (CCSDS) and it has been chosen for some deep 
space applications [9]. In the simulation, we set the noise 
variance of the contaminating heavy-tailed pdf (a 2) by 4.0 
and set the random walk signal of AWGN (Aa) by 10% of 
the noise variance of AWGN part (o~\). 

In this section, BER of the joint iterative channel estimation 
and decoding or the proposed approach is presented in order to 
compare with two methods. Conventional technique performs 



noise variance of AWGN part (af) estimation by second 
order moment SNR estimation based on hard-decisions of the 
previous codeword which they are feed back to multiply with 
that codeword to eliminate modulated information. Complex 
technique implements iterative SNR estimation which it give 
a high evaluation performance to estimate a\. 

Figure 6. depicts that the proposed technique has only 
about 0.2 dB at BER = 10~ 3 different from the SNR know 
and complex method. The efficiency of the proposed approach 
is close to the performance of ideal and complex techniques 
because the noise variance (o\ ) is re-calculated in every time 
of iterative LDPC decoder. Therefore, the noise variance (a\) 
has always been changed and obtains more accurate value in 
each time of LDPC iteration. In contrast, conventional method 
is the worst in performance compared with others. Even though 
having exceptional capability, the complex technique requires 
a large of processing time for that result depending on the 
number of SNR iteration. 



Bit Error Rate 




Figure 6. Bit-Error-Rate performance in the jiggling AWGN 
channel (the AR4JA (8192, 4096) code). 



TABLE I. makes comparison of the single block informa- 
tion for SNR estimation of those three techniques. Generally, 
look up table is used to find the answer of exponential 
operation. Hence, three kind of operations are taken into our 
consideration namely addition, multiplication, and look up 
table. 

Let length of single block information which pulsed RFI is 
erased is M. I is the number of iterative time of SNR estimation 
in the complex method, k represents the number of LDPC 
iteration. In the single block information, the conventional 
method requires the smallest operational numbers whereas the 
complex one demands a large amount of operation depending 
on number of I. The proposed approach has a little more 
computational complexity than the conventional technique 
depending on the number of LDPC iteration (k). Moreover, 
the complex method must perform look up table process while 
the others are absent in this process. 

For instance, in the multiplication operation, if the code 
length (N) and pulsed RFI area are equal to 10240 and 4% of 
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TABLE I. Comparisons of single block information for SNR 
estimation of three method. 





Conventional 


Complex 


Proposed 




method 


method 


method 


Addition 


M-l 


M+(2+3M)I 


(M-l)k 


Multiplication 


M+l 


3(1+I)+2M(1+2I) 


(M+l)k 


Look up Table 




4MI 





code length respectively, eliminated pulsed RFI signal (M) is 
equal to 10240-( 10240x0.04) « 9830. In the simulation, we set 
k and I by 50 and 100 respectively. According to TABLE L, 
the proposed approach succeed in decreasing the complexity 
in the multiplication operation by 87.56 % compared with the 
complex method. 

VI. Conclusion 

Joint iterative channel estimation and decoding has been 
proposed in this paper in order to mitigate the impact of 
pulsed RFI. The proposed approach presents joint iteration 
between channel evaluation and decoding technique so as to 
obtain more accurate the noise variance. The results show that 
the proposed technique has ability to attain less complexity 
computation than the complex method and enhance in the 
efficiency compared with the conventional method. Especially 
in the time varying SNR mismatch over AWGN channel. 
Therefore, the proposed approach is attractive to apply to 
various applications which pulsed RFI is involved with their 
requirements. 
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